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A novel theory Is developed for concentration polarization of non-Interacting particles in crossflow-filtration 
systems. This theory reveals that the extent of concentration polarization, as well as the behaviour of the per- 
meate flux, are characterized by an important dimensionless filtration number (W F = 4xajAP/3frT). There Is a 
critical value of the filtration number for a given suspension and operational conditions. When the filtration ' 
number is smal ler than the critical value, a polarization layer exists directly over the membrane surface and the 
wall particle concentration is determined by the pressure and temperature. At higher filtration numbers, a cake 
layer of retained particles forms between the polarization layer and the membrane surface. Mathematical 
models are constructed for both cases and analytical solutions for the permeate flux are derived. An Increase In 
permeate flux with increasing pressure Is predicted for all operational conditions. 



1. Introduction 

Crossflow filtration refers to a pressure-driven separation 
process in which the permeate flow is perpendicular to the 
feed flow. Separation in crossflow filtration is usually 
achieved by a membrane displaying different permeability to 
the solvent and solute (or particles). Crossflow filtration pro- 
cesses are conventionally classified as reverse osmosis 
(hyperfitaation), ultrafiltration, and inicrofiltration, depend- 
ing on membrane pore size and operation conditions. Cross- 
flow filtration has marked advantages over dead-end 
filtration and has been widely used to separate solutes from 
solutions or colloidal and particulate materials from suspen- 
sions. _ 
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Increases the , filter resistance and consequently reduces the v 
permeate flux through tl* niembra polar- 
ization is of considerable interest since a high permeate rate is 
most desirable in filtration processesTW a result, much 
research effort has been expended in this area. 

During the past two decades, numerous theoretical 
attempts have been made in an effort to understand concen- 
tration polarization and to predict permeate rate in 
crossflow-filtration systems. A representative achievement of 
these studies was development of the gel-layer model 1 * 2 This 
model has been used as a tool for directly interpreting experi- 
mental data in ultrafiltration and as a basis for further theo- 
retical development in microfiltration. In the latter case, 
shear-induced-diffuston and inertial lift are incorporated into 
the theory in order to obtain a better fit with experimental 
data. 3-5 Another popular theory in this field is the osmotic 
pressure model This model, however, has been shown to be 
fundamentally equivalent to the gel-layer model 6 

The gel-layer theory is derived from a single mass balance 
equation 1-3 with two unknowns: the permeate velocity and 
the particle concentration distribution over the membrane. 
Therefore, it is impossible to evaluate the permeate velocity 
from this theory without additional assumptions. The theory 
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becomes a semi-empirical model as it assumes a fixed surface 
particle concentration and adapts a mass transfer coefficient 
from theories of convective heat transfer to impermeable sur- 
faces. Neither has been theoretically proven valid for cross- 
flow filtration. Besides fundamental shortcomings, practical 
applications of this model are indeed limited. Although the 
gel-layer model has been developed to predict the so-called 
'tanking flux\ the theory cannot specify explicitly the criteria 
necessary for achievement of limiting flux. Furthermore, the 
model cannot be used for conditions where the permeate 
velocity is pressure dependent (&& before limiting flux is 
attained). 

The osmotic-pressure model has similar conceptual prob- 
lems. In this model, the wall particle concentration must be 
estimated from certain empirical relationships or simplified 
assumptions. This is one of the reasons why empirical expres- 
sions are often used to predict permeate flux in reverse 
osmosis systems. 7 ' 0 Evidently, this model cannot be applied 
to niicrofiltration and most ultrafiltration systems since 
osmotic pressure is negligible in these cases. 

In this paper, a novel theory for concentration polarization 
is developed, based on the hydrodynamics and thermodyna- 
mics of particle suspensions. It is shown that another funda- 
mental relationship (force or energy balance) is required, in 
addition to mass balance, to describe completely the concen- 
tration polarization in crossflow filtration. Similar to the gel- 
layer model, uniform non-interacting Chard 1 ) spherical 
particles are assumed in the mathematical derivation of the 
theory. It is expected that this new theory will provide a solid 
foundation for understanding and predicting permeate flux in 
crossflow-filtration systems. 



2. Physics in Crossflow Filtration 

Crossflow-filtration systems can have different configurations 
channel or tube). In the following sections, without loss 
of generality, a rectangular channel will be used in the theo- 
retical analysis and a complete rejection of the solute or par- 
tides by the membrane is assumed. Several fundamental 
physical principles are involved in the separation process by 
crossflow filtration. For instance, the flow field and drag force 
in crossflow filtration are described by basic theories of 
hydrodynamics and, more importantly, many bulk properties 
of particle suspensions in crossflow filtration are governed by 
thermodynamic principles. In this section, the important fea- 
tures of the flow field and driving force in crossflow filtration, 
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relevant to our theory, «c briefly w^jS^Jl* 
discussion on the thermodynamks of erossflow-filtration 

systems. 

XI FWwFIdd 

Flow processes over the membrane surface constitute^ major 
mechanism for concentration polarization m crossflow fil- 
tration. The transverse flow brings particles to the membrane 
surface and keeps the accumulated rjartieles m trie wanitt of 
the membrane surface. On in* other hand, the longjtudjnal 
(tangent) flow sweeps accumulated particles out of the filter 
and tends to reduce particle concentration on the membrane 

^Tteflow field in a rectangular channel with pamcable 
walKs) is governed, in r*inriple. by the Navier^tokes equa- 
tion, with the assumption of constant permeation velocity 
along the channel, the flow field in a channd wifc permeable 
walls was first obtained by Berman. 9 It shown that, in a 
thin layer near die permeable walls, the flow fiW at any loca- 
tion m the chaiind can be wcUar^ 
of a shear tangential flow and a constant perrjendicularfiow. 

As the suspension advances through the channel, the bulk 
fluid moves toward the inembrane surface owing to the 
applied pressure. The magnitude of this transverse flow is 
determined by the ar^ 

and the extent of concentration polarization. The P* 11 *^ 
velocity in crossflow filtration varies along tfiediaiuidlena* 
since the concentration-polarization layer thickness over the 
membrane surface increases gradually along the channel 
However, at a given location, the flow field ma tolaytr 
over the inembraiie surface can stffl be reasonably icpresent- 
ed by combination of a shear tongUudinal flow and a con- 
stant transverse flow. 

£2 Drh!agFo<c<s«i Resistance 

The driving force for permeate flow in crossflow filtration is 
the pressure drop from the bulk fluid tothe iwmeatestdeof 
the membrane (AP>. Particles or solutes are driven toward 
the membrane surface by the resulting rjressure gradients a 
closed system, as the permeate passes rbrc^ the mmbrane, 
the pressure gradient in the feed side the mmbrane wiB 
evetitualry vanish because of t^ 
arid the filtration rjrocessv/fflte 
atkm occurs in dead-end filtration. Crossflow filtration rep- 
resents an open system where the pressure gradient is 
maintained by a continuous supply of feed to the inlet and 
drainage of concentrated suspension from the outlet As a 
result, a constant permeate flow is attained. 

Permeate rate depends on the total resistance of the 
crossflow-filtration system. The total resistance to permeate 
flow is comprised of two components: (i) membrane resist- 
ance and (ii) resistance of the retained particles. Membrane 
resistance is usually a constant value, depending on mem- 
brane physical and chemical properties. Therefore, the flux of 
a pure solvent through the membrane, in the absence of 
retained particles, is proportional to the applied pressure. 
Resistance of the retained particles is attributed to the fnc- 
tional drag force resulting from the passage of permeate 
through the dense layer of particles which are stationary m 
the transverse direction. This resistance is a fundi n of the 
extent of particle accumulation and structure of the retained 
particle layer, and theref re depends on filtration per- 
formance. When the resistance is dominated primarily by the 
retained particles, filter performance is similar to that f an 
'ideal filter* as described below. 
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The Heal filter displays a complete rejection of solute or 
particles and shows no resistance to the process KqtnU fa 
such a filter, the entire pressure drop occurs in the retained 
r*rrkkfaya,tlu»greaty 

of concentration rxriarization. The concept of the Med filter 
wfll help us to understand the fundamental mceharrisnw of 
concentration polarization in crossflow-filtration systems. 

23 Principle of MhnmmnGibbs Energy 

In a crossflow-filtration system, because pressure decreases 
from the bulk toward the membrane surface, particles on the 
membrane surface wffl have minimum Gibbs energy. Par- 
ticles wfll accunwlate on the membrane surface because a 
systm tends to miimnize 

of particles in the vicinity of the surface is hindered by the 
thermal (random) motion of particles. The thermal motion of 
particles can be measured by the osmotic (thermodynamic) 
energy, which is a function of rjartide concentration The gra- 
dient of the osmotic energy is from the membrane surface to 
the bulk fluid. . . . . ^ 

At steady state, the total energy of a retained particle on 
any section of the filter is equal in the transverse dircctoi 
and there is no net movement of retained parades m tins 
Erection. The system is then said to reach a state of ther- 
(oodynamic equilibrium. In this case, the concentration <tf 
particles on the membrane surface is uniquely determined by 
the pressure drop applied to the accumulated particle layer. 
This describes the tliermodynamic foundation of concentra- 
tion realization in crossflow filtration. 

Because of the finite size of particles, the particle concen- 
tration on the membrane surface reaches its maximum value 
(closest packing) at a certain pressure. A further increase in 
the applied pressure beyond this point cannot be counter, 
balanced the osmotic ^ 

maximum value. This situation represents a rKWMqumTjnum 
(thermodynamic) system. In this case, particles m the bulk 
suspension have higher Gibbs energy than those on the mem- 
taamTsmface. The retained particles wfll move toward the 
membrane surface and form a cake layer between the polar- 
ization layer and the membrane surface. The cake fay* |»s » 
porosity corresponding to the maximum packing of the 
retained particles. Once a cake layer fonns, the pressure ^ 
across the polarization layer will be maintained at a critical 
value corresmmding to the porosity of a snaahnum packing. 
rL exc^k^r^essure will be absorbed by the cake 
layer. This process is Ascribed as cal» fonnatkm and can be 
considered as a special case of concentration pta" 8 *"- 

Since the Irydrodynarnics in the cake fayerc^cr from those 
in the concentration polarization layer, where r^de con- 
centration is below maximum packing, theories of crossflow 
filtration wfll be developed for both cases m the following 
sections. A criterion for cake formation wffl also be given 
along the development of the theories. 

& Concentration Polarization 

As liquid passes through the membrane, particles accumulate 
in the vicinity of the membrane surface and form a thin layer 
in which the particles are stationary in the transverse direc- 
tion with respect to the membrane surface. This stationary 
layer of retained partides provides additional reastence to 
permeate flow. The resistance of the retained particle layer 
singly depends on the total number of parades in this layer 
atKiTn theh^itial distribution. Therefore, this secto n wfll 
start with quantifying the number and concentration dis- 
tribution of retained partides in the concentration polariza- 
tion layer. 
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11 PtrtfcfcCooeeiitrtfion 

Since only the accumulated particles provide additional 
resistance to fluid flow, we use C to denote the excess number 
concentration of particles in the polarization layer. The 
actual particle concentration is C 0 + C, where C 0 is the bulk 
(feed) number concentration of particles. The distribution of 
accumulated particles in the stationary layer is governed by 
the following equations: 



dC 

p(x)C + i> — = 0 
dy 



(3.1) 



Cdy = M. 



where x and y are the longitudinal and transverse coordi- 
nates, respectively, c(x) is the permeate velocity at x, D is the 
particle diffusion coefficient, and M p is the total number of 
excess particles accumulated over a unit membrane suriace 
area. The direction of the permeate velocity here is defined as 
from the bulk to the membrane, so that i(x) is always posi- 
tive. The transverse coordinate y points into the flow from 
the membrane surface, to which we assign coordinate posi- 
tion y « 0l [In common formulations, the signs of the diffu- 
sive and convective terms in can. (3. 1 ) are opposite.] 

The solution of can. (3.1) with the mass balance condition 
eqn. (34) is 



C = *^cxr*-i<x)y/D] 
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The total number of retained particles, M p , is determined as 
follows. 

32 FirticHlmL Conservation 

At steady state, the particle flux along the channel is a con- 
stant and is equal to the particle flux entering (he inlet 
section of the channel. Therefore, one has the following flux 
conservation relationship: 



f 



ifCdv + tfCo0 = $>C o B 



(3.4) 



where u is the longittiainal velocity (here particles are 
assumed to take the velocity of the undisturbed flow), Uq and 
u are the average longitudinal velocities at the inlet section 
(jc«0) and at section of position x, respectively, and B is the 
channel height 

Inserting eqn. (33) into eqn. (3.4) and solving the resulting 
equation for M p> the total number of accumulated particles 
(per unit area) at position x, yields 



oCxOdY 



(3.5) 



where y is the shear rate of the longitudinal flow. Note that 
the shear longitudinal flow 

u = yy (3.6) 

and the fluid mass balance 



<«o-C)B= Wdx- 



were used to reach the above solution for Af 



(3.7) 



The particle concentrati n distribution in the channel is 
obtained by substituting eqn. (33) into eqn. (3 J): 



0*7 



[*x) J J^V) dv]exp[-c(x)y/D] (3.8) 
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It seems from eqn. (3J) that the wall particle conce ntra tion 
(O. at y m 0) is strongly dependent on the permeate velocity. 
However, as will be shown later, there is a more fundamental 
dependence of the wall concentration on thermodynamic 
variables. Even more striking is that the expression in the 
first square brackets of eqn. (18) is invariant in an ideal filter. 
This means that the wall concentration is not affected by the 
permeate velocity. We shall discuss this aspect later in more 
detail when the solution for the permeate velocity ts avail- 
able. 



(3-2) 33 Crossfl(nHoltrati(m Equation 



As mentioned before, there are two components for the resist- 
ance to permeate flow in a crossflow-filtration system, namely 
the membrane resistance and the accumulated particle layer 
resistance. The pressure drop through the accumulated parti- 
cle layer, AF P , can be calculated from 10 * 1 1 



FCdy 



(35) 



where F is the drag force exerted by a stationary retained 
particle on the permeate flow and £ p is the thickness of the 
accumulated particle layer. The pressure drop across the 
membrane is amply the product of the membrane resistance 
and the permeate velocity, gei^jhe^^^ 
c-rops across the membrane ^^accumulated particle 
laycr.ts equal to the applied pressure. This leads to the fol- 
lowing relationship: 

AF - Hx) + iM*)jc dF P.10) 

where AF is the pressure drop from the bulk to the permeate 
side of the membrane, A is the membrane resistance, ft is the 
Boltzmann constant, T is the absolute temperature and As is 
a c or rec tin g factor for Stokes* tew. The effective applied pies- 
sure drop, AF, is equal to the applied hydraulic pressure 
minus the osmotic pressure difference between the bulk sus- 
pension and the permeate. 

The expression in square brackets of eqn. (3.10) represents 
the drag force exerted by a stationary retained particle on the 
permeate flow. This expression is derived from Einstein- 
Stokes' tew, 11 combined with Happeft ceU model 11 to 
account for the effect of neighbouring retained particles. The 
term kTfD(=Gx(iaJ is the friction coefficient for a spherical 
particle of radius Op moving at tow speed in an unbounded 
fluid of viscosity p. The upper integration limit can be set to 
infinity because the concentration polarization layer is very 
thin compared to the channel height 

At high particle concentrations, as is often the case in the 
concentration polarization layer of crossflow filtration, the 
drag force exerted by a particle on the fluid deviates substan- 
tially from Stokes* law and a new relationship is required. 
This is usually furnished by introducing a correcting function, 
A s% into Stokes' expression. Happel's cell model 12 is used to 
determine in this study. According to the classical work of 
Happel 10 ' 12 



(3-11) 



where 0-{l- e)" 3 is a porosity dependent variable, with e 
being the porosity of the retained particles. 

The variable 0 is related to the number particle concentra- 
tion: 



fl«(iiwjC) 



(3.12) 



From eqn. (3*) and (3.12), one obtains 



90 



(3.13) 



Substituting cqn. (3.13) into cqn. (3.10) and changing the cor- 

' results in 



^ . m + £ J* j^yw* da (3.14) 

where 0 is the value of 0 at the membrane surface. Inserting 
cqn. (3.8) into cqn. (3. 12) and letting y « 0 yields 

^-[^■^r^r <3i5) 

Eqn. (3.14) is the governing equation for the permeate veloc- 
ity in crossflow filtration systems. 

34 Sofctftarf the Crosstab 

Applying the mean value theorem to eqn. (3.14) and substi- 
tuting cqn. (3.15) into the resulting equation gives 

AP - Mx) + ^ MP W*) 2 fa <**' C M> 



whercO<0*<0 w . 
Let us define 



0*7 



(3.17) 



Eqn. (3.16) can be then rewritten as 

AP « i«(x) + /I»|x|* £ sftf) <tf 

At x = 0. the integral in this equation vanishes, so that 

AP 

When x # 0, taking derivatives of both sides of eqn. (3.18) 
with respect to x results in 

(120) 



(3.18) 



(3.19) 



0 - £3 (1 + m*) fa + w** 3 



F!iipi"«Hn B the integration terms from eqn. (3.18) and (330) 
yields 

dx ~ 2AP-ai<x) 



(321) 
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filtration systems. By definition, the average flux, K, can be 
calculated from 



Eqn. (321) can be integrated with eqn. (3.19) as a boundary 
condition and the result is 

[(i) ,+ S x >» ,J - 3 s^- 4 =° (3J2) 

Eqn. (322) is a cubic algebraic equation of one unknown. The 
real root of the equation gives the general expression for the 
permeate velocity in crossflow-filtration systems: 

AP 

* x) ~(A 3 + 60AP a x) l ' J 

x <{[1W + 6/?AP J x) + 4]" 2 + 2} ,/J 
-|[Jl 1 i(jl 1 + 6fAP J x) + 4] WJ -2) w > (323) 

The average permeate flux is an important quantity which 
is commonly used to describe the behaviour of crossflow- 



if 



o(x)dx 



(324) 



where L is the channel kngth. A direct calculati on of ti ns 
integral is rather tedious because of the complex expression 
for *xl Fortunately, the value of this integral can be easily 
determined from eqn. (3.1 8). Thus, the average permeate flux 



is 



AP 



(325) 



35 Fatratiea Number NfandEsmmtlooofila^ 

In order to use eqn. (323) and (325). the value of AJMJnJl 
has first to be determined The estimation of 4^0*) can begm 
from eqa (3.14) by rewriting it as 

3*T AfV JL i-ia + i* 5 -* 6 

The term on the left-hand side is a dimenaonless number. We 
suggest to name it the filtration number for crossflow fil- 
tration and denote it as N F . That is 



! ^2eap» 

3*T 9 



(327) 



The filtration number, Np, can be considered as the ratio 
of the energy needed to bring a partide from tte membrane 
surface to the bulk suspension to the thermal (d"*""** 
energy of the partide. A larger value of N F indicates that 

more energy is needed to bring a F^^^^E"? 
surface to the bulk suspension. As will be shown later, N F is 
of fundamental importance to crossflow filtration. 

For most cases where concentration polarization is appre- 
ciably developed, h is safe to neglect membrane resistance, 
because it is usually much smaller than the retained parUcte 
layer resistance. Urns, the value of i«0») can be catoilated 
b> a* following procedure. When operational conditions ofa 
crossflow^ltration system are given, the filtiatoon number, 
N F . b estimated from eqn. (327) by substituting AP for AP.. 
The value of 0 V is then detennined from by solving tnc 
inverse problem of eqn. (326). Finally, by ar^ying the mean 
value theorem to eqn. (326) again, AJW * obUuncd from 

An accurate value of A^) b expected from the above 
procedure as long as the resistance of the filtration qsteme 
dominated by the concentration polarization layer. When tte 
membrane resistance is comparable with, or greater than, the 
resistance of the concentration polarization layer, this pro- 
cedure only provides an estimate of U<n However. the fil- 
tration number b usually very small and A^) * 1 when this 

situation occurs. . ^~Anr* 

Values of AJP), as well as 0 W , obtained by the procedure 
described above are Ibted in Table 1 for ranging from 0 
to 15.00. For a monodisperse suspension of rigid, spherical 
partides, the maximum packing occurs at a ™ £ 

0J6. The corresponding N F and 0 W values are 15.00 and 0^6. 
Therefore, in aU situations where concentrate nPp^za^ 
occurs, the AJfT) value can be directly obtained from Tabk 
1. When N F is greater than 15, the crossflow filtration 
becomes a n n^uilibrium system and a cake Uyer cm 
retained particles will form on the membrane surface, me 
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Trifel Values oT0 w and AJP) for different filtration numbers 
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Am 



0000 

i x i<r e 
i x ur 5 
aoio 
aioo 

O500 
IjOOO 
2.000 
3J000 
4000 
5.000 
6.000 
7.000 
8,000 
9.000 

taooo 

UJ000 
12J000 
13.000 
14J00O 
11000 



aooo 

0.002 
0.021 

ai98 

0381 
0556 
a633 
a705 
0.743 
0.768 
a786 
O800 
0U811 
0821 
0829 
0836 
0842 
0847 
08S2 
OSS6 
O860 



UOOO 
1J002 
IJ025 
IJ293 

I. 802 
2J910 
1946 
5.719 
7J25 
&842 
1O301 

II. 719 
11106 
14.468 
11809 
17.132 
18.441 
19.7i6 
21.021 
21295 
21559 



theory developed in this section cannot be used for such a 
system. A theory for non-equOibriuin filtration systems will 
be developed in Section 4. Schematic descriptions of concen- 
tration polarization and cake formation in crossflow fil- 
tration are presented in Fig. 1. 

Ranges of the filtration number for different crossflow fil- 
tration systems and typical operating conditions are listed in 
Table 2. It can be seen from the table that N F for reverse 
osmosis is much smaller than the critical value lor cake for- 
mation, while for inicrofiltration N F is much greater than the 



concentration polarization layer 




<6) 



suspension 



permeate flux v[x) 



ococeritrabwpotanzaton layer 




membrane 



T 

permeate flux v{x) 



Fig. 1 Schematic description of concentration pohrization and 
cake formation over a membrane surface in crossflow filtration, (a) 
Qmmrtrafr" polarization: there is a concentration gradient in the 
polarization layer. (6) Cake formation: particles accumulate in a cake 
layer at maximum packing. 




AP/Pa 

F? s I Permeate flux as a taction of applied pressure and bulk 
(feed) particle cc^cenlration. (a) linear scale; (6) enlargement of the 
initial stage and (c) logarithmic scale. The curves correspond to the 
following bulk (feed) cxmcentratkms: I, 1 * 10"; 2, 3 x 10"; 3, 
I x t0»;O x 10" and 5,1 x l©** m"*. Other parameters used in 
calculations: /= 109 Pa s m'\ r = «» S = 1S x ,0 ~" 
t«=05maodr = 293.15K. 

critical value. Only in some range of operating conditions of 
ultrafiltration, N f is of the same order of the critical value. 

X6 Behaviowiif the P^nneate Flux 

Average permeate fluxes as a function of pressure and feed 
particle concentration, calculated from eqn. (3.25), are plotted 
in Fig. 2 on both linear and logarithrnic scales. It can be seen 
from Fig. 2(a), which presents a set of typical permeate- 
pressure curves for crossflow filtration, that there are two 
stages in the operation of crossflow filtration. The permeate 
flux changes significantly with pressure in the first stage, 



Table 2 Ranges of the filtration number for Afferent types of crossflow filtration systems 



filtration system 



pressure/Pa (pst) 



particle size/m 



reverse osmesis 

ultrafiltration 

microfiltration 



lOMx 10* (150-600) 
W'-IO* (15-150) 
<3xl0*(<5) 



3.60 x 10" 10 
10-M0" 
5 x lO^-lO -5 



0.W9-O20 
0.I-10 6 
5 x I0M x 10 10 



3994 

while in the second stage, the flux increases much more mod- 
erately, at a rate which decreases as ttep 
The dependence of permeate flux on pressure can be better 

observed in Rg. 2(c). I^mmm^^^^i 

roughly as the first power of the preswierThis ineany^ 
concentration polarization is not significant in this range of 
applied pressure and the resistance to permeate flow is domi- 
nated by membrane resistance. This feature can be seen 
directly from eqn. (3.23). When X 3 > 6pAP 2 x (simflar to fil- 
tration of a pure solvent), we have 
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V = 



AP 



(3.29) 



(331) 



That is, the permeate flux is proportional to the applied pres- 
sure and to the inverse of the membrane resistance, 
tin the second stage, the concentration polarization is fully 
developed and the niter res i s tan c e is mainly detennined by 
the resistance of the retained partick la>er. This behaviour is 
wdl approximated by the ideal filter. For an ideal filter {U. 
I -» 0). eqn. (323) and (124) can be reduced to 

^= ( |WV' ^ x-««{^^)" , ("0) 

The flux changes as one third power of the pressure [if the 
change of A<jfi*) with pressure is not considered]. A similar 
derjendence of the average permeate flux on pressure in ultra- 
filtration systems has been recently reported. 13 

Tne eflect of partide concentration in the feed suspension 
on the permeate flux can also be seen from Fig. 2. At the 
stage where the permeate flux is controlled by concentration 
polarization, the flux decreases as the feed partide concentra- 
tion increases. 

XI Concentration Distrf^^ 

Combining eqn. (3iK (3.17) and (3.18) and rearranging gives 

AP 



[^I-^c(x)jexp(-r(x)r/D] (332) 

Eqn. (332) describes the partide concentration in the entire 
channel. Of particular interest is the parttcte concentration on 
the membrane surface (wall), at y « 0, Inspection of eqn. 
(332) reveals that the wall parttde concentration increases 
with x because the local permeate velocity i<jc) decreases with 
x. Iso-concentration curves in a crossflow filtration channel 
are presented in Fig. 3, on both linear and logarithmic scales. 
This figure indicates that the concentration polarization layer 
is very thin in the entire channd (of the order of 0.1 mm). 

When the permeate flux is mainly controlled by the resist- 
ance of the concentration polarization layer (i.e. 
1 J « 6/IAP 1 x) v the wall partide concentration is given by 



AP 



WW 



(3.33) 



Eqn. (333) states that when the membrane resistance can be 
neglected, the wall partide concentration of a given suspen- 
sion is determined solely by thermodynamic variables of the 
system (U. the pressure drop AP and the absolute tem- 
perature T). The wall particle concentrati n does not change 
with location and, even more striking, it is independent of the 
bulk partide concentration. Based on eqn. (333), it is worth- 
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Fm. 3 iso-coaceotritioo curves expressed as CfC m in a rectangular 
rh» nn M • ft fmear gale: (b\ Logarithmic scale. The aims uHicHwnd 
to the foOowtng CfC m values: 1, 095; 2. 0.70; X OSO; 4,030; 5,0.10 
and 6» 001. Parameters used in calculations: A« Mr Pa s m, 
7=100 s _, *C 0 «10** m"* ^,*15x MT 9 m, L = 0t5 m, and 
To 29X15 K. 

while to note that the expression in the last brackets of eqn. 
(330) and (331) is actually the ratio of the wall partick con- 
centration to bulk particle concentration. 

In the geMayer model, the wall particle concentration is 
assumed to be a constant independent of pressure for a given 
type of partkfes. It is obvious from the above analysts that 
this ftty«mpri<" «nnot he used, at least for the cases of pure 
concentration polarization, as otherwise it will violate the 
priKtplcsoftheraodynamic& 

4 Cake Formation 

When the filtration number, rV F , is higher than co. 15, for a 
monodtsperse suspension of rigid spherical particles, a cake 
layer of retained partides win form between the concentra- 
tion polarization layer and the membrane surface. By 'cake' 
we mean that the retained particles are packed so as to attain 
minimum porosity. In this section, a modd for the permeate 
flux will be developed for cases where the additional resist- 
ance of this cake layer is considered. This theory wfll com- 
plement the theory devdoped in the previous section, and it 
is particularly useful for most microfiltration and some ultra- 
filtration processes. 

41 Pressure Drop and Particle Distribution in the 
Polarization Layer 

When a cake layer forms, the entire pressure drop in a 
crossflow-filtration system occurs at the membrane, concen- 
tration polarization layer and cake layer, 



AP«AP„ + AP P + AP ( 



(4.1) 
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where AP a , AP, V and AP e are the pressure drops across the 
membrane, the polarization layer, and the cake layer, respec- 
tively. Since there is a cake layer under the polarization layer, 
the pressure drop through the concentration polarization 
layer remains at a critical value, which can be easily deter* 
mined by the following equation: 
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N rc kT 

i*4 



(4.2) 



where N Fc is the critical value of the filtration number at 
which the cake layer begins to appear. N R is a function of 
the closest-packing porosity. When the dosest- packing poro- 
sity of the retained particles is 0l36, N Fc is equal to 15.00. 

The distribution of retained parades in the polarization 
layer is similar to that without the cake layer [U. eqn. 
except for adjustment for the cake layer thickness in the y 



axis: 



^T^expt-BtxKy-^Z)] 



(43) 



where 6 9 aady-6 c are the thickness of cake layer and the 
distance from the cake surface at x, respectively. The particle 
concentration in the cake layer, on the other hand, is a con- 
stant which corresponds to the maximum packing of the 
retained particles. ' 

42 Cake Layer TOckness 

Applying Stokes" law and HappeTs cdl modd, the number of 
accumulated particles in the cake layer at position x, M e , can 
be obtained following a procedure similar to that described 
in the section on concentration polarization: 



OAP c 



(4.4) 



where 0^ is the 0 value for maximum packing. The cake 
layer thickness, S c , can then be calculated by 



0i 



(43) 



The expression in the square brackets of eqn. (45) is another 
dimensionless number. This dimensionkss number, denoted 
as N e , is termed the cake thickness fector. If. 



N ^A^OU 



(4.6) 



The physical meaning of this number win be dear from the 
following derivation, 

43 Pirtide-ftax Consecration 

Partidc-flux conservation along the filter channd at steady 
state yidds the following relationship: 



£"«Cdy = C 0 £ 



4*1 oV 



(4.7) 



The integral on the left-hand side of eqn. (47) represents the 
longitudinal flux of the retained partides in both the cake 
and polarization layers. 

An accurate calculation of this integral is difficult at 
present, because very little is known about the movement of 
retained partides in the cake layer and the effect of the cake 
layer on the fl w field above it It is assumed in this paper 
that the value of the left-hand side of eqn. (4.7) can be esti- 
mated through the following simplifying assumptions: (i) the 
vdodty of the accumulated particles in the cake layer is zero 
(immobile cake) and (u) the fluid field in the channd is undis- 



turbed by the cake layer. These assumptions comply with the 
intuitive anticipation that the build-up of the cake layer will 
enhance the downstream nwven^toftheitttmcdptrtkles. 

Based on the above assumptions, the integral on the left- 
hand side of eqn. (47) can be calculated and the eqiiation can 
be rewritten as 

Combining eqn. (4.5) and (48) to eliminate *5 C , and solving 
the resulting equation for M p yidds 



4x) 4xloV 



(49) 



9 Dyd+NJ 

It can be seen from eqn. (49) that the cake thickness factor, 
N c , reflects the effect of the cake layer thickness on the accu- 
mulation of retained particles in the polarization layer. 

44 RehtkMsUp between Pwneatt Velocity «d Cake 
Thickness 

Following a procedure similar to that used to derive eqn. 
(3.16), the permeate vdodty can be related to the pressure 
drop across the polarization layer by 

^^^[^^ (410) 

where 0 < 6** < 6L» - The value of M^tJ ean be found in 
Table 1. 

The solution of eqn. (410) gives the relationship between 
the local permeate vdodty 4x) and the cake layer thickness 
(represented by NJ: 

(411) 



or 



Furthermore, the amount of retained partides in the cake 
layer can be related to the permeate vdodty by combining 
eqn. (44H46) and eqn. (412): 



45 CrosshVm-fihiation Equation 

Considering the pressure drop across the membrane and the 
cake layer leads to the following relationship: 

AP _ AP p = AM + p£ 4KJ4x}]tf « (414) 

Substituting eqn. (413) into eqn. (414) results in 

AF* = i4*) + |P**4*) 3 ( 415 > 

where 
and 



16) 



D*7 



(4.17) 
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AP* can be taken, nominally, as the total pressure drop 
icross the membrane and the cake layer. 

46 Solultoacf theCrossflcw-fih^ 

The real root of cqn. (415) gives the general expression for 
the permeate velocity in a crossflow-filtration system where a 
cake layer forms: 



/AP*\" 3 



x <{IX WAJ* 2 *) + IT* + 1}"* 
- (PW^*) + IT* - U ,/3 > 

If the membrane resistance is negligible {U. AP e « AP 
- Afy, JV C can be calculated directly from eqn. (46). The 
local permeate velocity is then given by eon. (4.1U and the 
average velocity is determined by: 

(4.19) 

As N c -»a eqn. (411) and (4.19) are reduced to eqa (330) 
and (131). However, when N c > 1, (Le. N F > Npjt eqn. (4.11) 
and (4.19) become 

and 

Eqa(4Jfl) awl (431) cWer from eqa 
there are no pressure-dependent parameters in the former. As 
a result, the predicted permeate velocity from cqn. (4.13) and 
(114) wul increase strictly as one third power of the pressure. 

Trie detaMcxrjression for the tr^ 
can be obtained by combining eqn. (4.4), (45) and (418). 
When the membrane resistance can be n eg lected, the local 
cake layer thickness in the channel b simply given by 

Inspection of this equation reveals that as N F -»N*, j c -»Q. 
On the other hand, when N F >N R , 4« increases as two 
thirds power of the pressure. 
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11 Compkte filtration Cwnre 

It has been known for some time that a fanning flux may 
exist in crossflow filtration which is independent of the 
applied pressure. At the same time, however, different depen- 
dences of permeate flux on pressure have also been report- 
ed. 1 - 3 No sound explanation for this discrepancy has been 
offered yet Based on this study, it appears that the reported 
different dependences of permeate flux on pressure are a 
reflection of the different portions of complete filtration 
curves. 

A set of complete filtration curves, calculated based on the 
theory presented in Sections 3 and 4 is presented in Fig. 4. 
There are three regions in the curves, corresponding to the 
principal source of niter resistance. The first is a membrane- 
reststance dominated region, where the flux increases pro- 
portionally to the increase in pressure. This is followed by a 
polarization layer-resistattce dominated region, where the 
retained particles provide most of the resistance to permeate 
flow. In this region, the flux increases initially as one third 
power of the pressure, and the rate of increase de c lin es in the 
final stage because of the significant increase of AJ0*). The 
final region is a cake-layer-resistance dominated region, 
where the flux increases strictly as a cube root of the pressure. 

It is shown that there is no limiting flux in crossflow fil- 
tration when the assumptions used for the cake-layer analysis 
hold. Flux always increases with increasing pressure, but at 
different rates in different regions as shown in Fig. 4 The 
changes in ta beha^ 

c omp onen t of tte total ^resistance^jperrneate 

tte membrane to the cono 
(^tbxntoitejC&fi layer? In the last stage of the concentra- 
tion polarization resistance dominated region, the marked 
increase in A* counterbalances, to a certain degree, the effect 
of increase in pressure. Consequently, the dependence of flux 
on pressure in that stage is much smaller than one third 
power. This is possibly the reason ibrimrodudng the concept 
of hunting flux in previous studies of crossflow nTtratkm.'" 3 

S2 Happeft Ceil Mole! 

The hydrodynamics of concentrated suspensions are quite 
different from those of dilute suspensions, where the drag 
force exerted by a stationary spherical particle on moving 
fluid can be accurately calculated from Stokes* law. In a con- 
centrated particle suspension, HappeTs sphere-hvceU model 
is reported to predict the drag force with a considerable accu- 
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racy. 10 * 11 The vaMty of Happen model for electrolyte solu- 
tions (as for reverse osmosis) has sot, as yet, been reported fa 
the hteraturc. 

In constructing the filtration equation [U. eqn. (3.10)1 the 
key step involves relating the pressure drop across the polar- 
ization layer to the permeate velocity. This was a c co mplishe d 
by using Stokes* law combined with HappeTs model. This 
analysis, however, can be carried out in a general form as 
follows. 

Let be the chemical potential (per unit volume) of a com- 
ponent at concentration C. That is, 



*=/(Q 



(5.1) 



where C fa a function of the distance from the wall 
(membrane surface) due to concentration polarization. The 
osmotic (tbemodyiiamic) force per unit volume ts given b) 
the gradient of the chemical potential: 



dy 



<*2> 



When the component holds a stationary position with respect 
to the wall surface in the transverse direction, the osmotic 
force must be balanced by the applied pressure. Therefore, 
one has 



JL dCd* y 



M\ c 

Kit-*' 



where C w is the wall concentration of the coinponent and C* 
is a value between zero and C,. In the last step of the above 
derivation, eqn. (Xg) is used far C by letting y « a Eqn. (5J) 
holds for any chemical species Mate or particle) wfakh is 
retained by the membrane. 

Dilate suspensions of solute (or very small particles) are a 
special case where the osmotic energy per unit volume is 
given by ^ « kTC Therefore, we have 



d* 
dC 



(5.4) 



For t concentrated solution, a linear relationship does not 
exist between the osmotic energy art 
this case, the value of d^/dC can be determined by an 
approach tfmtiar to that used for partide suspensions, pro- 
vided the relation between chemical potential and solute con- 
centration [Le. eqn. (5.1H is known. For concentrated 
colloidal and particulate suspensions, Einstein-Stokes* law 
and HappeTs model are adequate to relate the chemical 
potential to the particle concentration, which results in 
d*/dC = ^0*)fcT. 

53 Units of PactkfeCoocestratioa 

Particle concentration is an important parameter in crossflow 
filtration. The volume fraction of particles is conventionally 
used to represent particle concentration in crossflow filtration 
but, as our theoretical analysis and ensuing discussion show, 
volume fraction is not a proper measure for particle suspen- 
sions. Most of the thermodynamic properties of suspensions, 
and hence the development of the concentration polarization 
layer in crossflow filtration, are related t the number con- 
centration f particles. At given operating conditi ns, suspen- 
sions with similar volume fraction but different particle size 
can result in quite different permeate velocities, and these dif- 
ferences cannot be explained simply by changes in the diffu- 
sion coefficient with particle size. The theory presented in this 
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paper dearly shows that the resistance of the cc^cc ntratfo n 
polarization layer is directly related to the number concentra- 
tion of retained particles through Stokes 9 tow combined with 
HappeTs model. Furthermore, with the number particle con- 
centration, some important conclusions can be directly 
derived from thermoo^namic principles. For example, when 
the wall concentration is below the saturation value, the 
number concentration, not the volume concentration, of par- 
ticks is uniqudy de t er min ed by the thennodynamtc proper- 
ties of the system. For the above reasons, the particle number 
concentration has been used in our theoretical development 

54 WaB-particfcCoMceatrttion 

The wall-partide concentration describes the extent of con- 
centration polarization and can be used to predict the per- 
meate vdodty in crossflow filtration. However, based on 
existing concentration-polarization theories, it is not dear as 
yet how to determine the waO-partide concentration. In the 
gd-layer theory, although the wall concentration is presumed 
to be a constant, no theoretical methods are provided to 
determine the constant value. On the contrary, the osmotic 
pressure theory takes the wall concentration as a variable, 
which changes with the applied pressure. However, an 
expression for the wall-particle concentration in the latter 
case is also not available. 

In this paper, we provide a dear and sound theory for cal- 
culating the wall-partide concentration. It has been shown 
that the waD-parucle concentration is untquery controlled by 
thermodynamic properties of the nitration system, c har acte r - 
ized by the filtration number N F . There is a critical value of 
the filtration number, N*, for a given suspension. When the 
N F of a crossflow-filtration system is greater than the of 
the suspension, a cake layer of fixed particle concentration 
(closest packing) will form on the membrane surface. In this 
case, the thickness of the cake layer, rather than the concen- 
tration in the cake layer, will change with the applied pres- 
sure. On the other hand, the wall-partide concentration 
changes with the appfied pressure when N v < N*. In this 
case, the wall-particle concentration is determined by the 
pressure drop appfied to the polarization layer. If the mem- 
brnoc resistance to permeate flow can be neglected, the wall- 
partide concentration is specified by eqn. (333). 

The possible existence of different wall-partide concentra- 
tions is a natural result of the theory presented in this paper. 
This theory can explain the mechanisms for the resulting 
wall-partide concentration in all crossflow filtration systems 
based on fundamental physical principles. It is, therefore, 
expected that this theory can be used to predict the behav- 
iour of crossflow filtration with more confidence. 

6. Sennmary 

A comprehensive theory for concentration polarization of 
non-interacting partides in crossflow filtration has been pre- 
sented in this paper. The essence of this theory is summarized 
below. 

The filtration number, AT F , defined by eqn. (3.27), is one of 
the most important parameters in crossflow filtration. When 
the closest packing of a given suspension is known, there is a 
critical value for N F . A cake layer will form between the con- 
centration polarization layer and the membrane surface when 
N F is greater than the critical value, while nly concentrati n 
polarization exists when N F is smaller than that value. Differ- 
ent sets of equations are formulated to describe the processes 
involved in these two cases. 

The filtrati n number has a simple dependence on particle 
size and applied pressure. For convenience, a working chart 
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10' W 
particle diameter/m 
fig. 5 A predominance diagram for concentration polarization and 
cake-fonnatkra regions as a function of applied pressure and particle 
* diameter 

is presented in Fig. S, where the regions of pure concentration 
polarization and cake formation are displayed as a function 
of pressure and partide diameter. A dosest packing porosity 
of 0.4 for uniform, spherical particles is assumed in construe- 
ting this diagram. 

From this chart, one can also find the appropriate equa- 
tions to use. In the pure concentration polarization region, 
the concentration distribution of accumulated particles in the 
channel is described by eqn. 032). The wall partide concen- 
tration is given by eqn. (333) when the membrane resistance 
is negligible. The local and average permeate fluxes are deter- 
mined by eqn. (3.23) and (3.25), and approximated by eqn. 
(3J0) and (3J1). respectively. 

In the cake-forming region, the local permeate flux is accu- 
rately given by eqn. (4.18). When the membrane resistance is 
negligible, the thickness of the cake layer is calculated from 
eqn. (422), and the local and average permeate fluxes are 
approximated by eqn. (4.11) and (4.19), respectivdy. When 
the cake layer resistance dominates, the focal and average 
permeate fluxes can be predicted by eqn. (420) and (421). 
respectivdy. 
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Foundation (Research Grant BCS-93081 19) and the State of 
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58385). The findings reported in this paper do not necessarily 
reflect the views of these agencies and no official endorsement 
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Notation 

As correction function for Stokes* law based on HappePs 

cell modd defined by eqn. (3.1 1) 
Op partide radius 
B channd height 

C excess partide number concentration 
C 0 bulk (feed) partide number concentration 
C w wall partide number concentration (at y = 0) 
O partide diffusion coefficient ;D « kTfintta P 
F drag f roe exerted by a stati nary particle on the per- 
meate flow 



ft 
L 
M t 

Mp 

Nc 

Nf 

Nfc 

P 

T 

ft 

u 

«0 
!<.X> 

V 

x 
y 

» 

f 
p 

/ 
/* 

e 
0 

0. 
0* 



AP 



AP e 

AP. 
Ar\ 



AP* 



I. CHEM. SOC FARADAY TRANS, 1995, VOL 91 

Boltzmarm constant, 13803 x 10"" J K" 1 
length of filter channd 

total number of partides (per unit area) accumulated 
m the cake layer at position x 
total number of partides (per unit area) accumulated 
in the polarization layer at position x 
cake-forming (actor defined by eqn. (4.6) 
dimensionless filtration number defined by eqn. (327) 
critical filtration number for cake forming 
pressure 

absolute temperature 

axial (longitudinal) fluid vdodty in the polarization 
layer; u = yy 

average tangential (cross) vdodty of fluid 
average inlet vdodty of fluid (at x = 0) 
local permeate vdodty (flux) 
average permeate vdodty (flux) 
axial (longitudinal) coordinate 
transverse coordinate 

fluid shear rate 

parameter defined by eqn. (3.17) 
parameter defined by eqn. (4.17) 
membrane resistance 
fluid viscosity 

cake or partide layer porosity 

parameter defined by eqn. (3.12); 0 = (1 - e) ,/3 

0 value corresponding to maximum packing of the 

retained partides 

the value of 0 at the membrane surface [y « 0) 

value of 0 between zero and 0 W 

chemical potential (per unit volume) 

thickness of cake layer 

thickness of concentration polarization layer 

effective pressure drop (applied hydraulic pressure 

minus bulk osmotic pressure) 

pressure drop across the cake layer 

pressure drop across the membrane 

pressure drop across the concentration polarization 

layer 

a nominal pressure drop defined by eqn. (4.16) 
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